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ABSTRACT 



The Supernova Remnant (SNR) IC 443 is an intermediate-age remnant well 
known for its radio, optical, X-ray and gamma-ray energy emissions. In this 
Letter we study the gamma-ray emission above 100 MeV from IC 443 as ob- 
tained by the AGILE satellite. A distinct pattern of diffuse emission in the 
energy range 100 MeV-3 GeV is detected across the SNR with its promi- 
nent maximum (source "A") localized in the Northeastern shell with a flux 
F = (47± 10) • 10~^ photons cm-^s"^ above 100 MeV. This location is the site 
of the strongest shock interaction between the SNR blast wave and the dense 
circumstellar medium. Source "A" is not coincident with the TeV source located 
0.4 degree away and associated with a dense molecular cloud complex in the SNR 
central region. Prom our observations, and from the lack of detectable diffuse 
TeV emission from its Northeastern rim, we demonstrate that electrons cannot 
be the main emitters of gamma-rays in the range 0.1-10 GeV at the site of the 
strongest SNR shock. The intensity, spectral characteristics, and location of the 
most prominent gamma-ray emission together with the absence of co-spatial de- 
tectable TeV emission are consistent only with a hadronic model of cosmic-ray 
acceleration in the SNR. A high-density molecular cloud (cloud "E") provides 
a remarkable "target" for nucleonic interactions of accelerated hadrons: our re- 
sults show enhanced gamma-ray production near the molecular cloud/shocked 
shell interaction site. IC 443 provides the flrst unambiguous evidence of cosmic- 
ray acceleration by SNRs. 

Subject headings: gamma rays: general — supernovae: general — supernovae: 
individual (IC 443) — cosmic rays — ISM: supernova remnants 
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1. Introduction 

Galactic cosmic-rays (CRs) are believed to be accelerated above lO^'^-lO^^ eV energies 
by powerful supernovae in our Galaxy (e.g., Shklovskii 1953; Ginzburg & Syrovatskji 
1964; Cesarsky 1980; Blandford & Eichler 1987), and gamma-rays above 70 MeV are 
expected to provide the crucial signature of hadronic acceleration. Shock-accelerated 
CRs interacting with the gaseous surroundings of Supernova Remnants (SNRs) produce 
gamma-rays by nucleon-nucleon interactions and neutral pion decay. Indeed, EGRET 
observations of SNRs have provided several important gamma-ray/SNR associations 
(ISturner Ez Dermer 1995| Esposito et al. 1996). However, the complex morphology of SNRs 



and the EGRET angular resolution prevented a definite resolution of this issue (e.g., Torres 
et al. 2003). In recent years, TeV detections of SNRs provided additional and very promising 
elements (e.g., Aharonian 2004; Aharonian et al. 2007a, 2007b, 2008; Yamazaki et al. 2006; 
Enomoto et al. 2006; Berezkho & Voelk 2006; Albert et al. 2007, 2008; Acciari et al. 2009). 
However, the hadronic interpretation of these detections usually requires a knowledge of 
SNR physical parameters or processes (e.g., electrons' Bremsstrahlung/inverse Compton 
emission vs. hadronic pion production, magnetic field strengths, the electron/proton 
number ratio for GeV-TeV kinetic energies, etc.). Difficulties in the interpretation of SNR 
gamma-ray and TeV data remain because of our poor knowledge of SNR distances and 
ages, ambient gas density, and problematic determinations of the non-thermal synchrotron 
and high-energy emissions. 

To unambiguously prove the CR acceleration by supernovae we need SNRs for which 
we can demonstrate that the ubiquitous and co-accelerated electrons do not dominate the 
observed gamma-ray and TeV emission. Given the variety of SNRs and the complexity of 
interactions with their environments, this task turned out to be very difficult to accomplishj^ 



^See also Butt 2009, for a recent review. 
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Current gamma-ray instrumentation can substantially improve this picture. Detailed 
gamma-ray mapping of SNRs can provide a first piece of evidence, i.e., confirming whether 
pion-generated emission in the energy range 100 MeV-10 GeV is mostly concentrated in 
sites where molecular clouds are strongly shocked by SNR expanding shells. As mentioned 
above, a second piece of evidence is required, showing that hadrons and not electrons 
are the main contributors to the detected gamma-ray emission. It turns out that IC 443 
provides both pieces of evidence. 



2. The Supernova Remnant IC 443 

The intermediate age SNR IC 443 (r ~ 10-20 thousand years) is located near the 
Galactic anticenter (1= 189.1, b = -1-3.0), and is at a relatively close distance from the 
Earth (~ 1.5 kpc) ( IWelsh fc Sallmen 2003p . It is one of the best studied SNRs because 
of its complex structure and interaction with its gaseous surroundings in the absence 
of strong diffuse Galactic emission. Radio flMufson et al. 19861 |Lehay 2004D , optical 



( IFesen fc Kirshner 1980P and X-ray (e.g., Petre et al. 1988; Asaoka & Aschenbach 1994; 
Kawasaki et al. 2002; Troja et al. 2006) mapping of the SNR show an asymmetric shape. 
The North-East rim expands in a relatively dense environment with a shock velocity 
of f s ~ 65 — lOOkms"^ (for an average number density of the unperturbed medium 
TT-i ~ 10 — 100 cm~^) ( IFesen fc Kirshner 1980p . and the South- West shell expands in a much 



more dilute medium (n2 ~ 1 — 10 cm~^). In the middle of the SNR, a dense and very massive 
(~ 1000 M0) molecular cloud complex in a form of a torus or "ring" surrounds from the 
exterior the expanding SN shell (Dickman et al. 1992, Lee et al. 2008). Carbon monoxide 
(CO) mapping of IC 443 shows several other smaller molecular clouds that interact with the 
expanding shell. In particular, the molecular cloud complex named "E" with a projected 
size of ~ Ipc and mass estimate equals to 23 Mq is the only prominent mass clump located 
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just in front of the expanding SN shell in the Northeastern rim (Dickman et al. 1992). 
Furthermore, detailed mapping of the J = 1-0 line of the formyl ion (HCO+, that traces 
compressed and heated gas within the SNR) clearly shows that both the massive "ring" 
and cloud "E" are physically interacting with the SN ejecta. Given the morphology of the 
SNR and its molecular clouds environment, IC 443 is therefore an ideal system to test the 
hypothesis of hadronic CR acceleration in SNRs. It provides a system which has all the 
required ingredients: a powerful SN of total explosion energy W ~ 10^^ ergs, both dilute 
and dense ambient gaseous environments surrounding the SN blast wave, and a complex of 
molecular clouds physically interacting with the SN shock. 



AGILE Gamma-Ray Observations of IC 443 



The AGILE satellite has been operating since 2007 April 23, ( ITavani et al. 2009|) 



The AGILE scientific instrumentation is very compact and is based on two co-aligned 
imaging detectors operating in the energy ranges 30 MeV - 30 GeV (GRID, Barbiellini 
et al. 2002, Prest et al., 2003) and 18-60 keV (Super- AGILE, Feroci et al. 2007). An 
anticoincidence system (IPerotti et al. 20061) . a calorimeter (ILabanti et al. 20^6|) . and a 
data handling system (Argan et al. 2004) complete the instrument. AGILE'S performance 



is characterized by very large fields of view (2.5 and 1 sr for the gamma-ray and hard X-ray 
bands, respectively), optimal angular resolution and good sensitivity (see Tavani et al. 2009 
for details about the mission and main instrument performance). 

During the first two years of operations AGILE observed several times the Galactic 
anticenter region. A total observing time of approximately 1 Msec was accumulated by 
AGILE, on IC 443, and a high-resolution mapping and spectral analysis in the energy 
range 100 MeV - 20 GeV were obtained using using standard AGILE gamma-ray selection 
procedures. Positional astrometry has been carefully checked by comparing the IC 443 
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Fig. 1. — Gamma-ray intensity maps (J2000 R.A. and DEC coordinates) above 400 MeV of IC443 
obtained by integrating all AGILE data. {Left panel:) gamma-ray intensity map above 400 MeV 
centered on IC 443. The color bar scale is in units of photons cm^^ s^^ pixel^^. Pixel size is 0.1 
degrees, and we used a 3-bin Gaussian smoothing. White contour levels of the gamma-ray intensity 
start from 0.0002 and increase in steps of 0.0000173. We also mark the positions and circular 
approximations of the 95% error boxes of the EGRET source 3EG J0617+2238 (green circle) and 
Fermi-LAT OFGL J0617. 4+2234 (black circle). The position of the TeV source associated with 
IC 443 is marked by a yellow circle and ellipse that give the 95% confidence level error boxes 
determined by MAGIC (|Albert et al. 2007D and VERITAS (|Acciari et al. 20091) . respectively (see 
also Butt 2009). The X-ray source CXOU J061705. 3+222127 is marked by a cyan circle. {Right 
panel:) optical image of IC 443 (Palomar Digitized Sky Survey) superimposed with the AGILE 
gamma-ray intensity contours above 400 MeV (same as left panel). The position of the TeV source 
is marked by the yellow circle (MAGIC) and ellipse (VERITAS). The angular distance between the 
centroids of the gamma-ray dominant source in the Northeastern rim and the TeV source is ~ 0.4 
degrees. The TeV source is located in the central part of the IC 443 at the location coincident with 
a massive molecular cloud complex (IDickman et al. 1992p . At the estimated distance of 1.5 kpc, 1 
arcmin subtends a distance of 0.44 pc. 



- 8 - 



gamma-ray emission with the nearby Crab and Geminga pulsars. The region is not 
substantially affected by the weak diffuse gamma-ray Galactic emission that is properly 
taken into account in our analysis ( IGiuliani et al. 2004^ . 



Fig. 1 shows the result of the AGILE detailed gamma-ray mapping (above 400 MeV) 
superimposed to an optical map of the Northeastern rim of IC 443. Diffuse gamma-ray 
emission with significant enhancements is detected across IC 443 in a pattern that closely 
resembles the SNR outer shell configuration. A most prominent gamma-ray enhancement 
(that we label "source A" in Fig. 1) is clearly detected in the Northeastern region of IC 443 
at a location coincident with the most active expanding SN blast wave. We also detect 
several other gamma-ray enhancements that follow a pattern surrounding the outer regions 
of the SNR shell. We label the remaining enhancements "B", "C", and "D" (see Fig. 1). 
The gamma-ray pattern of emission from IC 443 does not show any sign of significant 
variability. In this paper, we concentrate our analysis on the Northeastern rim of IC 443 and 
on source "A", postponing a more detailed study of the SNR to a forthcoming publication. 

Source "A" is detected at 11 sigma level at the Galactic coordinate location (l,b): 
(189.08,3.28) ±0.17 (stat.) ± 0.1 (syst.), with a flux F = (47 ± 10) ■ 10"^ photons cm'^ g-i 
above 100 MeV. We notice that the position of source "A" is consistent with that of the 
only prominent molecular cloud in the Northeastern region, i.e., cloud "E" of Dickman et 
al. 2002 which is in physical interaction with the expanding shell. Clearly marked in Fig. 1 
is the relatively small error box of the most prominent TeV source detected by MAGIC 
(lAlbert et al. 20071) and VERITAS (lAcciari et al. 2009^ . The pulsar wind nebula CXOU 
J061705.3-222127 (IC^aensler et al. 20061) is not associated with any prominent gamma-ray 



emission. 



Within the statistical and systematic uncertainties, the source "A" location is consistent 
with the EGRET (3EG J10617+2238, Hartman et al. 1999), AGILE (AGL J0617+2236, 
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Pittori et al. 2009) and Fermi-LAT (OFGL J0617.4+2234, Abdo et al. 2009) sources. 
While the EGRET source appeared to be marginally compatible with the position of the 
TeV source, the refined AGILE location is inconsistent at more than the 99 % confidence 
level with that of the TeV source. The location of the most intense TeV emission is indeed 
concentrated in the central part of IC 443 in apparent superposition with the massive 
molecular cloud "ring". The angular distance between source "A" and the TeV source 
centroid is more than 0.4 degrees, i.e., significantly larger than the location accuracies 
of AGILE (0.1 degree for this integration) and of MAGIC and VERITAS. We conclude 
that the positional difference between the gamma-ray source "A" and the TeV source is 
significant, and refiects the difference in the physical locations of the dominant 0.1-10 
GeV and TeV emissions, respectively. We infer a very important fact for the physical 
interpretation of our observations: the most prominent 100 MeV-10 GeV emission is not 
co-spatial with the TeV source. Furthermore, AGILE detects weak and diffuse gamma-ray 
emission from the central location of IC 443 in coincidence with the TeV source and the 
centroid of the massive molecular "ring". 

4. Discussion 

Let us denote by Xe,p the ratio between the electron and hadron number density 
normalizations (at momentum]? = 1 GcV/c). Electrons can emit high-energy photons by 
Bremsstrahlung and inverse Compton (IC) scattering. Emission by these two processes is 
unavoidably linked, in the sense that if an electron Bremsstrahlung contribution emerges as 
a prominent spectral component, a strong IC contribution is predicted in the TeV range. 
The same electrons emitting Bremsstrahlung gamma-rays (possibly enhanced by a dense 
cloud) also scatter unavoidably the ambient optical/IR photon bath of the SNR/interstellar 
medium and the cosmic microwave background (CMB) photons. A general conclusion 
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can be deduced for a typical SNR: a ratio Xe,p ~ 1 implies co- spatial 100 MeV and TeV 
emissions in SNRs. A detailed modelling of electron (and proton) emissivities of IC 443 
indeed predicts that for an ideal matter density of no ~ 1 cm~^ the IC contribution to the 
spectral power {vFi) dominates by a factor of ~ 10 over the Bremsstrahlung contribution 
for photon energies between 100 MeV and 1 GeV (Gaisser, Protheroe & Stanev 1998, 
GPS98). For higher energies, the IC dominance is even larger, leading to a predicted IC 
power at 300 GeV larger that a factor of ~ 100 than the Bremsstrahlung spectral power. If 
we apply these predictions to the specific case of the IC 443 Northeastern rim characterized 
by an external medium density ni ~ 10 — 100 cm~^, we expect the Bremsstrahlung power 
{^Fu)b near 0.1-1 GeV and the IC power {i'Fu)ic near 0.3-1 TeV to be approximately 
equal. This implies that a model with Xe,p ~ 1 predicts the electron-driven TeV emission 
to be co-spatial with the Northeastern rim, and in particular to be coincident with the 
0.1-1 GeV emission. Since this is not observed, we conclude that the SN shock must be 
characterized by a ratio Xe,p substantially less than unity. We note that a value Xe,p ~ 0.01 
has been indirectly deduced for IC 443 (GPS98) and other SNRs (e.g., Aharonian 2004), 
and in general agrees with what observed in direct CR measurements. 

The overall pattern of gamma-ray emission from IC 443, combined with the 
absence of detectable co-spatial TeV emission both near source "A" and over the 
whole semispheric shocked Northeastern rim, has far reaching consequences. Fig. 2 
shows the results of two leading models of gamma-ray emissivity produced by protons 
and electrons near source "A". We follow the treatment of Drury et al. 1994, and 
Berezhko & Volk 1997 for the normalization of the expected gamma-ray emission from 
accelerated hadronic CRs in terms of parameters of the late SNR evolutionary phase giving 
E^dN^/dE^ ~ 5 X 10"^^ (TeVcm-^s-^) {EsNfi-O^^ erg) (d/1 kpc)-^ m /, where EgN is the 
total supernova explosion energy, d the source distance, and / a geometrical factor. For 
simplicity, we assume that hadrons and electrons are accelerated at the Northeastern shock 
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with the same power-law energy index, i.e., N{E) = kE~". Models without energy cutoffs 
below 10 TeV for both electrons (-Ece) and protons (-Ec,p) contradict the observations. We 
consider then two model^ for source "A" as constrained by the data with Esn = 10^^ erg, 
d = 1.5 kpc, and / ~ 0.01. Model 1 is characterized by Xe,p = 0.03, a = 2.7, rii = 10^ cm~'^, 
Ec,e = 10 TeV, Ee,p = 200 GeV. Model 2 has Xe,p = 0.01, a = 2.1, rii = 10^ cm'^, 
i?c,e = 100 GeV, Ec^p = 100 GeV. Fig. 2 clearly shows that the hadron-produced gamma-ray 
emissivity has to be suppressed between 10 and 100 GeV. The electron Bremsstrahlung and 
IC contributions are also shown; both have to be consistent with the absence of co-spatial 
gamma-ray and TeV emission in the Northeastern rim satisfying the observational constraint 
z/F^(100MeV - lOGeV) > 100 uF^{0.2 - ITeV). We note that a value of Xe,p ~ 0.01 is 
consistent with the radio synchrotron emission observed in the Northeastern part of the 
IC 443 shell (Erickson & Mahoney 1985), for an average local magnetic field B ~ 10^^ G. 



5. Conclusions 

We obtain a satisfactory physical picture of IC 443 with very important implications 
for the problem of the cosmic-ray production by SNRs. Diffuse gamma-ray emission of 
hadronic origin is remarkably distributed in coincidence with the IC 443 outer shock. A 
clear gamma-ray enhancement is detected (source "A") in the Northeastern rim and close 
to the molecular cloud "E". This cloud provides an ideal target for nucleonic interactions 
of accelerated hadrons (protons and ions) with subsequent gamma-ray emission by neutral 
pion decay. Cloud "E" has a substantial target matter density for efficient hadronic 

^ In order to display the dynamic range of the most important parameter influencing the 
shape of the spectrum, we consider the values a = 2.1 (a realistic case), and a = 2.7 (a case 
with no spatial diffusion). 



-VI- 




Fig. 2. — The AGILE gamma-ray spectrum (between 100 MeV and 1 GeV) and the TeV upper 
hmits of the region centered at the source "A" of IC 443 compared with two different models for the 
gamma-ray emissivities. Solid curves: hadronic production of gamma-rays by neutral pion decay 
(blue curve: model 1, black curve: model 2). Dashed curves: Bremsstrahlung emission by electrons 
(blue curve: model 1, black curve: model 2). Short- dash- dot curves: electron IC emission on the 
cosmic microwave background (blue curve: model 1, black curve: model 2). The MAGIC upper 
limits are assumed to have values 3 times smaller than the reported flux levels of the TeV source 
at the center of IC 443 (jAlbert et al. 20071 lAcciari et al. 20091) . The AGILE spectral data agree 
with those reported by EGRET ( |Esposito et al. 1996[ ). 
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interaction with a timescale Tpp ^^^o ~ (5 ■ 10^ yrs) (n/10^ cm~^)^^. This timescale turns out 

to be comparable with the age of the remnant, and imphes a high efficiency of gamma-ray 
production by the SN blast wave-cloud interaction. By assuming an average particle kinetic 
energy of order of 0.1—1 TeV, our gamma-ray observations of IC 443 and the proposed 
physical interpretation are in agreement with the total energy in cosmic-ray particles being 
1-10% of the total estimated SN explosion energy of about 10^^ ergs. 

At the source "A" site, both the proton and electron distributions have an effective 
energy cutoff of order of 100 GeV. We note that several models of particle acceleration in 
high-density regions close to the acceleration site consider the possibility of effective energy 
cutoffs of order 0.1-1 TeV (e.g., Drury et al. 1996, Baring et al. 1999, Malkov et al. 2002). 
These effective cutoffs can be caused by a combination of inefficient acceleration at larger 
energies and particle diffusion. The CR maximum kinetic energy Em produced by shock 
acceleration can be estimated as D{Em) = RgVs/k', where D{E) is the particle spatial 
diffusion coefficient, k' ~ 30 for the late evolutionary phase, and Rg and Vg are the shock 
radius and speed, respectively (e.g., Berezhko 1996). An upper limit for Em applicable for 
source "A" of IC 443 can be obtained with a Bohm approximation for D, and -R^ = 1 pc, 
Vg ~ 10* cm~^ s~^, and i? ~ 5 x 10^^ G. We obtain Em ~ 1 TeV, indicating that the active 
SNR shocked region currently contains only relatively low kinetic energy particles with the 
highest energy CRs diffused away a long ago. Our observations of the intermediate-age 
IC 443 support this important aspect of the particle acceleration mechanism. 

It is also interesting to address the lack of prominent gamma-ray emission in coincidence 
with the TeV source in the high-density medium at the center of the SNR. A possible 
explanation is based on a combination of energy dependent particle diffusion and nucleonic 
interactions in the very dense central medium reached by the SN shock (ITorres et al. 2008|) . 
In any case, the absence of detectable TeV emission from most of the SNR argues for a 
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sub-dominant electron contribution to the emission. 

We presented here two specific emission models in agreement with the high-energy 
observations and yet predicting different fluxes in the TeV range. Future deep TeV 
observations of IC 443 may reveal a weak TeV emission associated with the Northeastern 
rim and the source "A" region. Fig. 2 provides an example of how these observations may 
help in further constraining the theoretical models. 

IC 443 turns out to be the first SNR clearly providing evidence for hadronic cosmic-ray 
acceleration and interaction with its gaseous surroundings. A leptonic model of emission is 
in contradiction with the combined gamma-ray /TeV observations of the Northeastern rim 
of IC 443 and in particular of the source "A" region. A hadronic model of emission agrees 
in a natural way with our observations, and confirms the hypothesis that supernova blast 
waves can efficiently accelerate protons and possibly other ions in our Galaxy. 

We thank an anonymous referee for his/her comments that led to an improvement of 
the discussion of our results. The AGILE Mission is funded by the Italian Space Agency 
(ASI) with scientific and programmatic participation by the Italian Institute of Astrophysics 
(INAF), and the Italian Institute of Nuclear Physics (INFN). This investigation was carried 
out with partial support from the ASI contract n. 1/089/06/2. 
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